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Abstract 
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We investigate the propagation of a charged particle in a spatially constant but time 
dependent pseudoscalar background. Physically this pseudoscalar background could be 
provided by a relic axion density. The background leads to an explicit breaking of Lorentz 
invariance; as a consequence processes such as p —¥ pj or e — > ej are possible within some 
kinematical constraints. The phenomenon is described by the QED lagrangian extended 
with a Chern-Simons term that contains a 4-vector which characterizes the breaking of 
Lorentz invariance induced by the time-dependent background. While the radiation induced 
(similar to the Cherenkov effect) is too small to influence the propagation of cosmic rays 
, in a significant way, the hypothetical detection of the photons radiated by high energy 

cosmic rays via this mechanism would provide an indirect way of verifying the cosmological 
relevance of axions. We discuss on the order of magnitude of the effect. 

> 

o 

oo 
in 

CO 

d 
o 

> 

• rH 

X 

S3 



ICCUB-10-062 
UB-ECM-PF-10-36 

*e-mail: domenec.espriu@cern.ch 
^e-mail: arencer@gmail.com 



1 



1 Axions 



Cold relic axions resulting from vacuum misalignment [TJ [2] in the early universe is a popular and 
so far viable candidate to dark matter. If we assume that cold axions are the only contributors 
to the matter density of the universe apart from ordinary baryonic matter its density must be [3] 

p ~ i(r 30 gcm- 3 ~ l(T 46 GeV 4 . (1) 

Of course dark matter is not uniformly distributed, its distribution traces that of visible matter 
(or rather the other way round). The galactic halo of dark matter (assumed to consist of axions) 
would correspond to a typical value for the density [1] 

p a ~ l(T 24 gcm- 3 ~ l(T 40 GeV 4 (2) 

extending over a distance of 30 to 100 kpc in a galaxy such as the Milky Way. Precise details 
of the density profile are not so important at this point. The axion background provides 
a very diffuse concentration of pseudoscalar particles interacting very weakly with photons 
and therefore indirectly with cosmic rays. What are the consequences of this diffuse axion 
background on high-energy cosmic ray propagation? Could this have an impact on cosmic ray 
propagation similar to the GZK cutoff [5]? This is the question we would like to address here. 

The fact that the axion is a pseudoscalar, being the pseudo Goldstone boson of the broken 
Peccei-Quinn symmetry [B], is quite relevant. Its coupling to photons will take place through 
the anomaly term; hence the coefficient is easily calculable once the axion model is known 

AC = g ail ^-jFF. (3) 

Two popular axion models are the DFSZ[7j and the KSVZ[8j ones . In both models <7 a77 ~ 1. 
Here a is the axion field and f a is the axion decay constant. Further details are provided in 
section 3. 



2 Cosmic Rays 

Cosmic rays consist of particles (such as electrons, protons, helium and other nuclei) reaching 
the Earth from outside. Primary cosmic rays are those produced at astrophysical sources (e.g. 
supernovae) , while secondary cosmic rays are particles produced by the interaction of primaries 
with interstellar gas. In this work, the effect of axions on the propagation of these cosmic rays 
will be studied. We will separately consider proton and electron cosmic rays and ignore heavier 
nuclei because the effect on them will be far less important as will become clear later (the 
axion- induced Bremsstrahlung depends on the mass of the charged particle) . 



2.1 Cosmic Ray Energy Spectrum 



We are interested in the number of protons in cosmic rays. Experimentally, one sees that the 
number of cosmic ray particles with a given energy depends on energy according to a power law 



J(E) = NiE~ 



(4) 



where the spectral index 7^ takes different values in different regions of the spectrum (see [9]). 
For protons we have 



J P (E) 



5.87 • 10 19 £- 2 - 68 10 9 < E < 4 • 10 15 

6.57 • 10 28 £~ 3 - 26 4 • 10 15 < E < 4 • 10 18 

2.23 • 10 16 £T 2 - 59 4 • 10 18 < E < 2.9 • 10 19 

4.22 • 10 49 E~ 4 - 3 E > 2.9 • 10 19 



(5) 
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while for electrons the power law isjlUj 

5.87 • 10 17 £T 2 - 68 E < 5 • 10 10 
4.16 • 10 21 £- 3 - 04 E > 5 • 10 10 



Je(E) 



(6) 



and the flux typically two orders of magnitude below that of protons, although it is more poorly 
known. Our ignorance on electron cosmic rays is quite regrettable as it has a substantial impact 
in our estimation of the radiation yield. 




Energy [eY] 



Figure 1: Proton cosmic ray energy spectrum 



Note that the above ones are values measured locally in the inner solar system. It is known 
that the intensity of cosmic rays increases with distance from the sun because the modulation 
due to the solar wind makes more difficult for them to reach us, particularly so for electrons. 
In addition, the hypothesis of homogeneity and isotropy holds for proton cosmic rays, but not 
necessarily for electron cosmic rays. Indeed because cosmic rays are deflected by magnetic 
fields they follow a nearly random trajectory within the Galaxy. We know that on average a 
hadronic cosmic ray spends about 10 7 years in the galaxy before escaping into intergalactic 
space. This ensures the uniformity of the flux, at least for protons of galactic origin. On the 
contrary, electron cosmic rays travel for approximately 1 kpc on average before being slowed 
down. However, because I ~ t 1 ^ 2 for a random walk, 1 kpc corresponds to a typical age of an 
electron cosmic ray ~ 10 5 vr|llj. In addition, the lifetime of an electron cosmic ray depends on 
the energy in the following way 

t(£)^5xl0 5 (^)yr = |, (7) 

with To ~ 2.4 x 10 40 . To complicate matters further, it has been argued that the local interstellar 
flux of electrons is not even representative of the Galaxy one and may reflect the electron debris 
from a nearby supernova ~ 10 4 years ago[12j. 



2.2 The GZK Cut-off 

The GZK (Greisen-Zatsepin-Kuzmin) limit [5 J states that the number of cosmic rays above a 
certain energy threshold should be very small. Cosmic rays particles interact with photons from 
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the Cosmic Microwave Background (CMB) to produce pions 

7CMB + P — >p + tt° or 7CMB + P — >n + ir + . (8) 

The energy threshold is about 10 20 eV. Because of the mean free path associated with these 
reactions, cosmic rays with energies above the threshold and traveling over distances larger than 
50 Mpc should not be observed on Earth. This is the reason of the rapid fall off of the proton 
cosmic ray spectrum above 10 20 eV as there are very few nearby sources capable of providing 
such tremendous energies. 

Note that the change in slope of the spectrum at around 10 18 eV is believed to be due to 
the appearance at that energy of extragalactic cosmic rays. 



3 Solving QED in a Cold Axion Background 

In this section we shall describe in great detail the theoretical tools needed to understand the 
interactions between the highly energetic cosmic rays we have just described and the cold axion 
background described in the first section. 

The interaction of axions and photons is described by the following piece in the lagrangian 

-£(177 = 5(177 y - -^ 1 F^vi (9) 

where 

fv» = ^ ua ^F aP (10) 

is the dual field strength tensor. 

The axion field is originally misaligned and in the process of relaxing to the equilibrium 
configuration coherent oscillations with q = are produced, provided that the reheating tem- 
perature after inflation is below the Peccei-Quinn transition scale [6]. In late times the axion 
field evolves according to 

a(t) = aocos(m a t), (11) 
where the amplitude ao is related to the initial misalignment angle. With this, Q becomes 

£ a77 = 5a 77 ;^--^-ao cos(m a t)F^ ' F^ v = g ail -^a cos(m a t)e^ al3 d^A v F al3 . (12) 

^ 7r la 7T J a 

Integrating by parts (dropping total derivatives) and taking into account that e^ ua ^ d^F^ = 0, 
we get 

am a a . ijk 
£a77 = ga 7 7 7 — sm(m a t)e J AiF jk , (13) 

7r J a 

where Latin indices run over the spatial components only. 

A cosmic ray particle (which travels at almost the speed of light) will see regions with quasi- 
constant values of the axion background, of a size depending on the axion mass, but always 
many orders of magnitude bigger than its wavelength. Thus, we can approximate the sine in 
(|13p by a constant (^, for example). Then, it can be written as 

£ a77 = \^A V F^\ (14) 
where ry^ = (r/, 0,0,0) and r/ = 4g ary The "constant" r\ changes sign with a period 

71 J a 

~ l/m a . 
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The oscillator has energy density p a = -dj^ ax = -{m a ao) 2 , so m a ao = \/2p a . Then, the 



constant rj is 



4a yfZp^ 20 
V = 5a 77 -, — ~ 10 eV, 



(15) 



V" /a 

for Pa = 10~ 4 eV 4 and f a = 10 7 GeV = 10 16 eV. 

The extra term in (|14[) corresponds to Maxwell-Chern-Simons Electrodynamics. Although 
in Maxwell-Chern-Simons Electrodynamics one can have in principle any four- vector 77^, the 
axion background provides a purely temporal vector. We shall assume rf" to be constant within 
a time interval l/m a . 



3.1 Euler-Lagrange Equations 

In the presence of an axion background the QED Lagrangian is 

C = --F^F^ + ^ (iff - e 4 - m e ) ^ + \miA^ + \^A V F^ 



(16) 



Here also an effective photon mass has been considered (equivalent to a refractive index, see 
[2]). It is of order 

(17) 



m 7 ~ Ana 



n. 



The electron density in the Universe is expected to be at most n e ~ 10 7 cm 



10 



-21 



:V 3 . 



This density corresponds to m 7 ~ 10~ 15 eV, but the more conservative limit (compatible with 
[T3] ) m 7 = 10~ 18 eV will be used here. 

The second term of (|16p gives the kinetic and mass term for the fermions and also their 
interaction with photons. Dropping it, we get the Lagrangian for (free) photons in the axion 
background (see [2] for further details): 



-\f^F^ + im^A" + \r)v.A v F» v 



1 



1 



d^A v {d^A v - d u A») + + - A e^^A v d a A p . 



(18) 



The Euler-Lagrange (E-L) equations are 



da 



dC 



dC 



d(d a A x ) dA } 



(19) 



dC 
dA, 
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K \m 2 y g^A^A u + - A e^^A u d a A p 

' m-irr' !.. 4- r;'" 1 I ,, 1 4- ^ 



ra:,(g Xu A u + 5 "%) + L^ Vfl d a Ap 



x 2 
2 ' t 



1 



m^ A + 



■q^daAp. 



(20) 



0a 
0a 
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d(d a A x 
i 



-(3 CT ,4 A - A A ff ) + -e^V^, 



-0 CT CT .4 A + A (T A' J + Je^r^A,. 



1 



(21) 
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Rearranging the indices, the equations are 



DA X + d x d a A a - m 2 A x - ^e^^daAp = 0. 



(22) 



If we choose the Lorenz gauge d a A a = the second term vanishes. The equations can also be 
written as 

- g^UAp - g xp m 2 A p - \e^ a ^d a A p = 0. (23) 

We are interested in writing these equations in momentum space. To this end, define the Fourier 
transform of the field: 



The relevant derivatives are 
and 

The E-L equations are then 

d A k 
(2vr) 4 



d a A/3 



UA, 



d 4 k 

d 4 k 



-ikx 



Mk). 



(2vr) 4 

d 4 k 

(2vr) 4 



{-ik a )e- ikx Ap{k) 
(-k 2 )e- ikx Ap(k). 



^A (fc 2_ m 2 )+ *^ ka 



e~ ikx Ap{k) = 0. 



Therefore, 



9 ^{k 2 -m 2 1 )+ l -e^ a r lll k a 



M k ) = o, 



or 



K^A v {k) = 0, = g^(k 2 - m 2 ) + -e^r, a kp. 

3.2 Polarization Vectors and Dispersion Relation 

We now define 

S\ = e^rj a k^ Xpa ri^. 



(24) 

(25) 
(26) 

(27) 

(28) 
(29) 



(30) 



This can be put in a more convenient form using the contraction of two Levi-Civita symbols 

,0123 v 



^^Xpa^ va ^ = — 3!(5^5p <5^j (the minus sign is there because in Minkowski space eoi23 = — e C 



S y.u = 



fa • k) 2 - 7] 2 k 2 } g^ -(vk) (ifk v + k^r] u ) + k 2 ifrf + r] 2 k»k u . 



It satisfies 



SW = S»X = 0, S = S» = 2 ( V .k) 2 - V 2 k 2 , s^s uX = -sl. 



If if = (r], 0, 0, 0) we have S = 2i] k > 0. Now we introduce two projectors: 



flak p. 



(31) 



(32) 



(33) 



These projectors have the following properties: 

p£> rh/ = p£> kv = 0t g^ = l, (P^)*=I^ = P^, 
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p^p ±Xu = p^ p^p tXu = o, p^+pr = -s^. 



(34) 



With these projectors, we can build a pair of polarization vectors to solve ()29|) . We start from 
a space-like unit vector, for example e = (0, 1, 1, l)/\/3. Then, we project it: 



= P£ u e u . 

In order to get a normalized vector, we need 



(35) 



g/2e% + ?? 2 (e- fc) 2 _ _1 (e • k) 2 
S ~~2 + 2 P 

(this is of course negative because e is space- like). Then, the polarization vectors are 



(36) 



(e-fc) s 



2/fc 2 



-1/2 



y 7 — e±4f 
These polarization vectors satisfy 

gfip£± £± — —I, d\iv£± ^zp — u 

and 

2 

4*4 + 44* = --r = — ^r- 



With the aid of the projectors, we can write the tensor in (|29p as 



= g^(k 2 - m 2 ) + a/| (PT - ) 



(37) 

(38) 
(39) 

(40) 



Then we have for = (w±, fc) 



(A: 2 - m 2 ) =F 



= (fc 2 -m 2 Tr7|£|)4 = ( <4 ~ ^ ~ "4 =F ^\k\ ) 4- (41) 



Therefore, = 4 is a solution of fljgj} iff 



u±(k) = \/m 2 ± rj\k\ + k 2 . 



(42) 



This is the new dispersion relation of photons in the cold axion background in the approximation 
where rj is assumed to be piecewise constant. 



4 The Process p — > p^ 
4.1 Kinematic Constraints 

We now consider p(p) — > p(q)'j(k), or e(p) — > e(q) , j(k). This process is forbidden in normal 
QED due to the conservation of energy. It is, however, possible in this background (the cold 
axion background even allows the process 7 — > e + e~ , see [IS])- Momentum conservation means 
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q = p — k. Calling m the mass of the charged particle (proton or electron), conservation of 
energy leads to 



E(q) + uj(k) = E(p), y m 2 + (p — k) 2 + y m 2 ± rj\k\ + k 2 = \J m 2 + p 2 , 
y/E 2 + k 2 - 2pkcos6+ ^m 2 ± r/k + k 2 — E = 

In the last line, a lighter notation has been adopted: 

E = E(p) = \J m 2 +P 2 , p = \p\, k = \k\, p-k = pkcos9. (44) 

As will be seen, if r, is positive (negative) the process is only possible for negative (positive) 
polarization. Therefore, ±i] = — \r>\ in these cases. To take into account both of them, we will 
use the minus sign and write r\ instead of \rf\. 
Squaring twice yields 

(4E 2 - 4p 2 cos 2 9 + Apr] cos 9 - r] 2 )k 2 - 2(2E 2 r] + 2m 2 p cos 9 - m 2 r])k + (4E 2 m 2 - m*) = 0. (45) 
Neglecting m 7 , rj in front of m, E this is: 

{E 2 - p 2 cos 2 9+pr] cos 9)k 2 - [E 2 r] + m 2 p cos 9)k + E 2 m 2 = 0. (46) 
This equation has two solutions 



E 2 r\ + pm 2 cos 9 ± EJE 2 r\ 2 — 4E 2 m 2 + 4p 2 m 2 cos 2 9 — 2pm 2 r] cos 9 



k± 2(E 2 -p 2 cos 2 9 + pr, cos 9) ' ^ 

These solutions only make sense if the discriminant A is positive. With the approximation 
cos 9 ~ 1 — \ sin 2 9, the condition A > is 

■ 2n^ [v 2 tf -2pm 2 r, + m 2 (j} 2 -4m 2 )\ 

sm 9 v^d - g • (48) 



Which can be rewritten as 

sin 2 # < 

4p 



2 ^ r, 2 1 



7 4p 

where 



S " r ° ~ 4v 2 m 2 1 - ^> ~ P+){p ~ (49) 



P± = ^ ± ^ Jl - c ± 2 -^lJl - J?L (50) 
7] 7] W 4t7i~ 7/ W 4m^ 



It is clear that p + > and p_ < 0. For sin 2 9 to be positive we need 



2/71777/ / 77^ 

P>P+=m = 1 \ 1 -T-2- ( 51 ) 

77 y 4m^ 

This is the threshold below which the process cannot take place kinematically. The energy 
threshold (E 2 h = m 2 + p 2 h ) is: 

E, h = (52) 
V 

When 77 — > 0, the threshold goes to infinity (as is expected: the process cannot happen if 77 
vanishes) . 
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There is another relevant scale in the problem: m 2 /rj. It is many orders of magnitude above 
the GZK cut-off. Therefore, we will always assume the limit p <C m 2 /r]. The maximum angle 
of emission for a given momentum is given by (|49|) : 

rj 2 1 

sin 2 6 max (p) = ^-(p-p+){p-P-)- (53) 

4^ 777^ 1 - 4^ 

Its greatest value is obtained when p is large (p Pth)'- 

sin 2 max = (54) 

Since this is a small number, photons are emitted in a narrow cone ^ max — — . This justifies 
the approximation made for cos#. 

At #max(p)) the square root in (|4"T|) vanishes and 

mwp)] = Me^GO] = Ptfe<<p<<m2/r? ) ^ (55) 

The minimum value for the angle is 9 = 0: 



E 12 ?? + pm 2 ± (^E 2 r] 



pm 2 



MO) ^ ^5 s (56) 

This gives the maximum and minimum values of the photon momentum. In the limit pth <C 
p -C m 2 /r] they are: 

^max = fe+(0) = (57) 

2 

k min = k_(0) = ^. (58) 
77 

These two values coincide at the energy threshold. 

Here we can see that the process is possible for negative (positive) polarization only if 77 > 
(rj < 0). Otherwise, the modulus of the photon momentum would be negative. 

Note that the incoming cosmic ray wavelength fits perfectly within the X/m a size, so it 
indeed sees an almost perfectly constant r\. Whether 77 is positive or negative there is always a 
state with slightly less energy to which decay and lose part of its energy (of 0(r/)) emitting a 
soft photon. So even if the process is a rare one it does not average to zero. An exact analysis 
will we presented elsewhere. 



4.2 Amplitude 

The next thing we need is to compute the matrix element for the process. Using the standard 
Feynman rules we get 

iM = u{q)ie^u(p)el{k). (59) 

Its square is 

\M\ 2 = u{q)ie^u{p)el{k)[u{q)ie 1 v u{p)el{k)}* = e 2 el{k)e v {k)tr [u{q)u{q)^u( P )u{p)Y] ■ 

(60) 

We now must sum and average over initial and final proton helicities, respectively. We do 
not average over photon polarizations because the process is possible only for one polarization. 
Performing the trace: 

W = \e 2 el{k)e p {k)tr[{4+m)^{p' + m)Y] 
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IxW 7 - 




Figure 2: The solution k± of the energy conservation equation (|43p . It can be seen that ([5 
and (1581) are indeed the maximum and minimum values. 



= -e 2 e;(k)e u (k)tr[MY + ^l u ] 

= \e%{k)e v {k)W - q a p a g^ + q^ + m 2 g^\. 

Using 4- momentum conservation, (|38p and the fact that p a p a = m 2 , we get 

pWF = 2e 2 [-p a k a + 2 £ * £ „pV1 = 2e 2 [ p a k a + + e M eJ) pV 

Now we use (f39|) to get 

The averaged square amplitude is then 



^ — = p sm i 

rj z k z 



\M\ 2 = 2e 2 (-p a k a + p 2 sin 2 0) . 



(61) 



(62) 



(63) 



(64) 



The first term is positive: 



m 2 — rjk — 2Euj 1 
—p a k a = — Eco+pk cos 9 = —Euj — pk — - 



1 



2pA; 



l r] (k--^) = ^(k-k min )>0, (65) 
2 T] 2 



so |.M| 2 is clearly positive. 

4.3 Differential Decay Width 

The differential decay width is 



d T = (27r) 4 6^(q + k-p)^\M\ 2 dQ, 



where the phase space element is 



dQ 



d 3 c 



d 3 k 



(2vr) 3 2£(g) (2vr) 3 2w(A;) ' 



(66) 



(67) 
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We can use 5^(q + k — p) to eliminate d?q. The remaining 8 is the conservation of energy. 
Therefore, E(q) = E — oj. Next we use a property of the Dirac delta function 



where x\ are the zeros of the function. In our case, we consider E(q) a function of cos#: 



(68) 



5[E(q)+u-E] = 6( y / E 2 + k 2 -2pkcos9 + uj-E 

-2pk ' ' 



2yjE 2 + k 2 - 2pk cos a 



5 cos 9 



m 2 — rjk — 2Euj 
-2pk 



(69) 



Next we write Sk = k 2 dkd(cos 9) dip, integrate the (p angle (factor of 2ir) and use the delta to 
eliminate d(cos9). This fixes the value of cos 9: 



cos 9 = 

Finally, the differential decay width is 

m a k 
dT 



2 Epu 



m 2 — rjk — 2Elo 
-2pk ' 



( — p a k a + p sin 2 9)dk, 



(70) 



(71) 



where a = e 2 /4-7r and sin 9 is given by (|7Up . This decay width can be written more conveniently 
for future computations: 



dT a 1 



dk 



k<jj 



[A{k) + B(k)E^ 1 + C{k)E~ 



,E 2 r\ 



k), 



(72) 



with 



A(k) = 4(r)k-m 2 ), B(k) = Auj{m 2 -r,k), C(k) = -2m 2 k 2 +2r]k 3 -m't-ri 2 k 2 +2m 2 r]k. (73) 



4.4 Effects on cosmic rays 

We now want to compute the energy loss of protons in this background 



dE _ dtdE _ 1 
dx dx dt v 



LodT 



(74) 



Using the previous results and v = p/E, the energy loss is (with the integration limits given by 
fl£3) and (I55D) 



dE 
dx 



a 1 

2p^ 
a 



kdk 



^min 



-(rjk — m 2 ) + p 2 (l — cos 2 9) 



dk 



h ■ 



2r]k 2 - (Am 2 + 2m 2 + r] 2 )k + 2r)(2E 2 + m, 



(4E 2 + m 2 )- + AErjJm 2 - rjk + k 2 + AEm 2 



a 



rj 3 E 6 

m 6 7? 3 



1 



m 2 — i]k + k 2 
k 



(4m 2 + 2m 2 + rf 



rfE*_ _ 
m 4 rj 2 
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+2 V (2E 2 +m 2 ) [5^-^] -m 2 (4£ 2 + m 2 )ln ( ) + 
The leading term is 

dE _ a 2r/ 2 £ 4 ar/ 2 ^ 2 ^ ar/ 2 ^ 2 
8p 2 m 2 4m 2 v 2 Am 2 

The energy as a function of the traveled distance is then 



(75) 



(76) 



E(x) = ^ . (77) 



The fractional energy loss for a cosmic ray with initial energy E(0) traveling a distance x is 

E(0) - E(x) %&E{0)x 



E(0) " 1 + ^E(0)x 



(78) 



4m 2 

2 



This loss is more important the more energetic the cosmic ray is. However, 4^ is a very small 
number. If we take .E'(O) = 10 20 eV (the energy of the most energetic cosmic rays) and x = 10 26 
cm (about the distance to Andromeda, the nearest galaxy, therefore larger than the galactic 
halo) the energy loss is smaller than 1 eV. For less energetic cosmic rays, the effect is even 
weaker. 

As we have seen, the effect of the axion background on cosmic rays is quite negligible. 
However, the emitted photons may be detectable. Using m 7 = 10~ 18 eV and r\ = 10 -20 eV 
as indicative values and having in mind the GZK cut-off for protons (and a similar one for 
electron^) the emitted photon momenta fall in the range 

1(T 16 eV < k < 100 eV (79) 

for primary protons and 

10~ 16 eV < k < 400 MeV (80) 

for primary electrons. 

The number of cosmic rays with a given energy crossing a surface element per unit time is 

d 3 N = J(E)dEdSdt , (81) 

where J(E) is the cosmic ray flux. These cosmic rays will radiate at a time t. The number of 
photons is given by 

d 5 N y = d 3 N dT( - E /^ dkdt = J(E) dT ( E,k ^ dEdkdt dSdt. (82) 
dk dk 

Assuming that the cosmic ray flux does not depend on time, we integrate over to obtaining a 
factor t(E): the age of the average cosmic ray with energy E. Since we do not care about the 
energy of the primary cosmic ray (only that of the photon matters), we integrate also over E, 
starting from E m i n (k), the minimum energy that the cosmic ray can have in order to produce 
a photon with momentum k, given by (|57|) . Therefore, the flux of photons is 



d 3 iV 7 f°° dT(E,k) 



' dEt{E)J{E) ^^ nj , Eth = 2 r ^. (83) 



dkdSdt JE min (fc) dk r] 



x It is very doubtful that electrons could be accelerated to such energies but it is irrelevant anyway for the 
present discussion as the intensity is extremely small at these energies 
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Next we assume that t(E) is approximately constant and take t(E) ~ T p = 10 7 yr for protons 
and t(E) ~ T e = 5 • 10 5 yr for electrons. We know that this last approximation is not correct as 
t(E) ~ 1/E but at this point we are just interested in getting an order of magnitude estimate 
of the effect. 

The photon energy flux is obtained by multiplying the photon flux (|83|) by the energy of a 
photon with momentum k: 



I(k) = uj(k) 

(k)>E t h 

aT r °° 
~8k 



dE t(E)J(E)^ 



dE Ni 



A{k)E~^ + B{k)E~^ +l) + C(k)E~ 



(7i+2) 



(84) 
(85) 



where E m i n (k) = m^J |, see (|57p . Numerically, the only relevant term in the decay rate is 4r/k, 
from A(k). The integral can then be approximated by 



T(U\ ar]T J i E min(k)] E min (k) ,-2=1 



In 



1 



The value "f m i n is to be read from ([5|) or © depending on the range where E m i n (k) falls. 

Substituting the numerical values we obtain the following approximate expressions for I p (k) 
and I e (k) 



I p (k) = 6 x 



I e (k) = 200 x 



Tn 



10 7 yr 



V 



10- 20 eV 
V 



1.84 



A: 



10- 7 eV 

2.02 / k 



-0.84 



m 2 s 1 sr 1 . 



-1.02 



m 2 s 1 sr 1 . 



(87) 



5 x 10 5 yry V10- 20 eV/ V 10 " 7 eV , 

As mentioned above these expressions are only indicative and assume constant average values 
for the age of a cosmic ray (either proton or electron). For a more detailed discussion we 
encourage the reader to examine our recent paper [16J . From this latter work we include the 
following figure describing the radiation yield 



Uk]( 



2-1 

m s sr 



') 



*p y + e-*e y 




k(eV) 



Figure 3: Radiation yield using the exact formulae and a more appropriate parametrization 
of the electron cosmic ray average lifetime as a function of the energy. From |16] , Note that 
electrons in general dominate the effect at low energies. 
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5 Conclusions and Outlook 



In this work, the effect on charged particles of a mildly (compared to the particle momentum) 
time dependent pseudoscalar background has been investigated. We have been interested both 
in proton and electron cosmic rays. 

This effect is calculable because the axion background induces a modification of QED that 
is exactly solvable. This modification has some interesting features, such as the possibility of 
the photon emission process p — > pj and e — > ej (which we have termed as axion-induced 
Bremsstrahlung processes). Kinematical constraints on the process have been reviewed, in 
particular it is seen that it is only possible for proton energies higher than a certain threshold. 
The energy loss of protons in such a background has been computed. For protons that survive 
the GZK cutoff this loss is totally negligible. 

However, the radiated photons could still be detected. Their flux and energy spectrum have 
been computed in some detail. Since the energy threshold depends on the mass of the charged 
particle, it is lower for lighter particles. Also, the energy loss is proportional to the mass squared 
of the charged particle, so the effect is more important for electrons. The value of k m [ n does not 
depend on the charged particle mass, so the radiated spectrum is no very different for electrons 
or protons (however the average lifetime of electron and protons cosmic rays is quite different 
and this has an observable effect on the power spectrum of the radiation). 

We refer the interested reader to |16j for a more comprehensive description of this phe- 
nomenon and on the possibility of this diffuse radiation being measured. We summarize however 
the main conclusions below. 

The dominant contribution to the radiation yield via this mechanism comes from electron 
(and positron) cosmic rays. If one assumes that the power spectrum of the cosmic rays is 

_ 7-1 

characterized by an exponent 7 then the produced radiation has an spectrum k 2 for proton 

primaries, which becomes k 2 for electron primaries. The dependence on the key parameter ry ~ 

■Jo 1+7 2+7 

—jr^ comes with the exponent r\ 2 and rj 2 for protons and electrons, respectively. However 
for the regions where the radiation yield is largest electrons amply dominate. We have assumed 
that the flux of electron cosmic rays is uniform throughout the Galaxy and thus identical to the 
one observed in our neighbourhood, but relaxing this hypothesis could provide an enhancement 
of the effect by a relatively large factor. The effect for the lowest wavelengths where the 
atmosphere is transparent and for values of r/ corresponding to the current experimental limit is 
of 0(1O -1 ) mJy This is at the limit of sensitivity of antenna arrays that are already currently 
being deployed and thus a possibility worth exploring. 

In the case of radiation originating from our galaxy the main unknown in the present discus- 
sion is whether the flux of electron cosmic rays measured in our neighbourhood is representative 
of the Galaxy or not. Since it is possible to relate this flux to the galactic synchrotron radi- 
ation one could deduce the former from measuring the latter. It appears [T7] that either the 
total number of electron cosmic rays is substantially larger than the one measured in the solar 
system, or the galactic magnetic fields have to stronger than expected. This issue remains to be 
further quantified. No attempt has been made to quantify the signal from possible extragalactic 
sources either. 

One should note that the effect discussed here is a collective one. This is at variance with 
the GZK effect alluded in the first section - the CMB radiation is not a coherent one over 
large scales. For instance, no similar effect exists for hot axions. A second observation is that 
some of the scales that play a role in the present discussion are somewhat non-intuitive (for 
instance the 'cross-over' scale rrii/r] or the threshold scale m 7 m p /r/). This is due to the non 
Lorentz-invariant nature of this effect. Finally, it may look surprising at first that an effect that 
has such a low probability may give a small but not ridiculously small contribution. The reason 
why this happens is that the number of cosmic rays is huge. It is known that they contribute 
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to the energy density of the Galaxy by an amount similar to the Galaxy's magnetic field[18j. 

There are several aspects of the present analysis that could be improved to make it more 
precise, particularly a piecewise constant oscillating axion background, or one with a serrated 
time profile for that matter, could be solved easily without having to appeal to special functions 
(the sinus profile involves Mathieu functions). This will be presented elsewhere but the present 
analysis suffices to indicate the order of magnitude of the effect. 

We hope that the present mechanism help to assess the presumed relevance of cold axions 
as a dark matter candidate. 
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